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Abstract 

Jim Propp's rotor router model is a deterministic analogue of a random 
walk on a graph. Instead of distributing chips randomly, each vertex serves 
its neighbors in a fixed order. 

Cooper and Spencer (Comb. Probab. Comput. (2006)) show a remark- 
able similarity of both models. If an (almost) arbitrary population of chips 
is placed on the vertices of a grid Z"* and does a simultaneous walk in the 
Propp model, then at all times and on each vertex, the number of chips 
on this vertex deviates from the expected number the random walk would 
have gotten there by at most a constant. This constant is independent 
of the starting configuration and the order in which each vertex serves its 
neighbors. 

This result raises the question if all graphs do have this property. With 
quite some effort, we are now able to answer this question negatively. 
For the graph being an infinite fc-ary tree (fc > 3), we show that for 
any deviation D there is an initial configuration of chips such that after 
running the Propp model for a certain time there is a vertex with at least 
D more chips than expected in the random walk model. However, to 
achieve a deviation of D it is necessary that at least exp{Q,{D^)) vertices 
contribute by being occupied by a number of chips not divisible by fc at a 
certain time. 



1 Introduction 

The rotor-router model is a simple deterministic process first introduced 
by Priezzhev et al. and later popularized by Jim Propp. It can be viewed 
as an attempt to derandomize random walks on graphs. So far, the "Propp 
machine" has been studied primarily on infinite grids Z**. There, each vertex 
a; G Z'^ is equipped with a "rotor" together with a cyclic permutation (called 
a "rotor sequence") of the 2d cardinal directions of Z''. While a chip (particle. 
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coin, . . . ) performing a random walk leaves a vertex in a random direction, in 
the Propp model it always goes in the direction the rotor is pointing. After a 
chip is sent, the rotor is rotated according to the fixed rotor sequence. This 
rule ensures that chips are distributed highly evenly among the neighbors of a 
vertex. 

The Propp machine has attracted considerable attention recently. It has 
been shown that it closely resembles a random walk in several respects. The 
first results were due to Levine and Peres Q (and later Landau and Levine ) 
who compared random walk and Propp machine in an aggregating model called 
Internal Diffusion-Limited Aggregation (IDLA). 

Cooper and Spencer [I[ compared both models in terms of the single vertex 
discrepancy. Apart from a technicality, they place arbitrary numbers of chips on 
the vertices. Then they run the Propp machine on this initial configuration for 
a certain number of rounds. A round consists of each chip (in arbitrary order) 
performing one move as directed by the Propp machine. For the resulting chip 
arrangement, they compare the number of chips at each vertex with the expected 
number of chips that the random walk would have put there, starting from the 
same initial configuration and running for the same time. Cooper and Spencer 
showed that for all grids Z"*, these differences can be bounded by a constant Cd 
independent of the initial setup (in particular, the total number of chips) and 
the run-time. For the case d = 1, that is, the graph being the infinite path, the 
optimal constant ci is approximately 2.29 0, for the two-dimensional grid it is 
C2 « 7.87 i. 

This raises the question whether the Propp machine on all graphs simulates 
the random walk that well. In this work, we show that the infinite fc-regular tree 
behaves different. We prove that here arbitrarily large discrepancies can result 
from suitable initial configurations. However, to obtain a discrepancy of D, at 
least exp(f2(Z?^)) vertices have to participate by being occupied by a number of 
chips not divisible by k at some time. 

While the work cited above and ours in this paper primarily aims at under- 
standing random walks and their deterministic counterparts from a foundations 
perspective, we would like to mention that meanwhile the rotor router mecha- 
nism also led to improvements in applications. An example is the quasirandom 
analogue of the randomized rumor spreading protocol to broadcast information 
in networks 

2 Preliminaries 

To bound the single vertex discrepancy between the Propp machine and a ran- 
dom walk on the fc-regular tree we first introduce several requisite definitions 
and notational conventions. 

Let G = (y, E) be the infinite fc-regular tree, also known as the Cayley tree 
and the Bethe lattice. We fix an arbitrary node to be its origin 0. |x| denotes 
the distance between the origin and vertex x. 

Wc first describe the Propp machine in detail. Each vertex x is equipped 
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with a rotor pointing to one of its neighbors. If A is the 'direction' the rotor is 
pointing to, then we denote this neighbor by x + A. This notation is in analogy 
to the one used in the grid case. Though we keep the additive notation, here 
the directions would rather be the generators of the nonabelian group generated 
by the expression \df = 1), whose Cay ley graph is an infinite fc-regular 

tree. Denote by DIR := {di, . . . , dk} the set of generators ("directions"). 

Each vertex also holds a rule describing how the rotor is moved. This rule is 
encoded by a cyclic permutation of the directions called rotor sequence. Again, 
if the rotor of vertex x is pointing in direction A, then we denote by NEXt(A) := 
NEXt(A,x) the new direction after one move. Since there will be no danger of 
confusion, we will always omit the explicit mention of the vertex. Nevertheless, 
we should stress that we do not require all vertices to have the same rotor 
sequence. 

Chips move according to the rotors. That is, if at a certain time t there is a 
single chip on vertex x (with rotor pointing to A) , then at time t + 1 this chips 
is on vertex x + A, the neighbor the rotor pointed to, and the rotor is updated 
to point in direction NEXt(A). If there are more than one chip on a vertex, 
they all move in one round, but applying the rotor principle one after the other 
(in an arbitrary order — we do not care about individual chips but only about 
the number of chips present on each vertex) . 

With these rules, the Propp machine describes a process that is fully deter- 
mined by the initial setting of the chips and the positions of the rotors. For all 
X G F and < G No, we denote by /(x, <) the number of chips on vertex x and 
ARR(x, t) the direction of the rotor associated with x after t steps of the Propp 
machine. Hence f{-,0) and ARR(-,0) describe the initial configuration, which 
determines all other values of / and ARR. 

For completeness, let us briefly review the random walk model. Here, we do 
not have rotors. Instead, in each round, each chip independently and uniformly 
at random chooses a neighbor of his current position and moves to there. This 
process is described by the initial setting of the chips together with all random 
decisions made. 

However, since we shall be only interested in the expected number of chips 
on the vertices, we may instead regard the following linear machine Here, 
in each time step every vertex splits its pile of chips evenly. When a pile of 
£ chips splits evenly on some vertex, £/k chips go to each neighbor. By the 
harmonic property of random walks, the (possibly non-integral) number of chips 
at vertex x at time t is exactly the expected number of chips in the random 
walk model. 

Note that the linear machine again is a deterministic process fully described 
by the initial numbers of chips on each vertex. Given such an initial configura- 
tion, we denote by i5(x, t) the (fractional) number of chips at time t on vertex 
X. Note that, again, £'(•,0) is just this initial configuration. Note further that 
£;(x, t) = i Eagdir + a, t - 1) for aU i e N and x e 

To compare the Propp machine with the linear machine, we start both with 
identical settings of the chips. Hence we have /(x, 0) = £'(x, 0) for all x G V. A 
configuration describes the current state of machine we regard. A configuration 



3 



of the Propp machine assigns to each vertex x e F its current (integral) num- 
ber of chips and the current direction of the rotor. For the hnear machine, a 
configuration is a simple mapping V — > Q>o, describing the (fractional) number 
of chips on each vertex. 

As pointed out in the introduction, there is one limitation without which 
neither the results of [iHl] nor our results hold. Note that since G is a bipartite 
graph, chips that start on even vertices never mix with those starting on odd 
vertices. It looks as if we are playing two noninteracting games at once. How- 
ever, this is not true. Chips in different bipartition classes may affect each other 
through the rotors. We therefore require the initial configuration to have chips 
only on one class of the bipartition. Without loss of generality, we consider only 
even initial configurations, i. e., chip configurations supported on vertices at an 
even distance from the origin. 

To analyze the behavior of the linear machine, we need the following no- 
tation. By H(x,t) we denote the probability that a chip from a vertex with 
distance x to the origin arrives at the origin after t random steps ( "at time i" ) 
in a simple random walk. Then, 

H{x,t) = k~*n{x,t) (1) 

with n{x,t) counting the number of paths of length t between two vertices 
at distance x on the infinite /c-rcgular tree. It is easy to verify the following 
properties of n{x, t): 



n(0,0)==l, 

n{x,0) ==0 for all a: > 1, 

n{0,t) =kn{l,t~l) foralH>l, 

n{x,t) =7i(x - 1,< - 1) + {k - l)n{x + l,t - 1) for aU x,t>l. 
Finally, we write x ~ t to mean that |x| = t (mod 2). 



3 Mod-fc-forcing Theorem 

For a deterministic process like the Propp machine, it is obvious that the initial 
configuration (that is, the location of each chip and the direction of each ro- 
tor), determines all subsequent configurations. The following theorem shows a 
partial converse, namely that (roughly speaking) wc may prescribe the number 
of chips modulo k on all vertices at all times by finding an appropriate initial 
configuration. An analogous result for the one-dimensional Propp machine has 
been shown in [2]. 

Theorem 1 (Mod-fc- forcing Theorem). For all initial directions of the rotors 
and any tt: x No — > {0, 1, . . . , (fc — 1)} with 7r(x, t) = for all x 9^ there 
is an initial even configuration /(x, 0) that results in subsequent configurations 
satisfying /(x, t) = 7r(x, t) (mod fc) for all x and t > 0. 
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Proof. Fix arbitrary initial directions for the rotors. Let us call /, or more 
precisely, an initial configuration (/(x, 0))xev fine for a subset S CV x No, 

• if it yields /(x, t) = 7r(x, t) (modfc) for all (x, G S, and 

• if it is even, that is, for all x G wc have /(x, 0) = if x 7^ 0. 

Consider first the initial configuration defined by /'°-'(x, 0) :— 7r(x, 0) for all 
X e y. Clearly, /(o) is fine for V x {0}. 

Assume now that there is a T > and an initial configuration f'^'^^ which 
is fine for V x {0, ...,T}. In order to obtain an /'■"^^^•' that is fine for 
V X {0, ...,T+ 1}, we shall modify /(^) by defining (x, 0) := /(^)(x,0) + 
e(x) k'^^^ for appropriately chosen e(x) e {0, 1, . . . , (fc — 1)} with e(x) = if 
X7^ 0. 

Observe that a pile of k'^'^^ chips splits evenly T + 1 times. Therefore, for 
all choices of the £(x), /(^+i)(x,t) ee (modfc) for all (x,t) with t < T. 

This implies that /C^+i) is fine for x {0, . . . , T}. 

We shall define e inductively such that Z*-"^^^^ is also fine for V x {T + 1}. 
Let e^°^(x) := for all x G 1^. Assume that for some > 0, e^^^ is such that the 
initial configuration /^'^^ (x, 0) := /^^^ (x, 0) (x) fc^+i is fine for all (x, T+l) 
with |x| < 9. Trivially, this is fulfilled for = 0. For each x with |x| = 9, we 
choose some y(x) with |y(x)| = T + 1 + 9 such that the distance between 
X and y(x) is T + 1. Note that this implies that the mapping x i-> y(x) is 
injective. We define e^^+^^ as follows. For all x with |x| = 9, let e'>^+^^ (y(x)) := 
(7r(x,r + 1) - /W(x,r+ 1)) modfc. For all z € ^ \ {y(x) : |x| = 9), we set 
£(f+i>(z) £^^^(z). We first argue that is an even initial configuration. 

As /^^^ is an even initial configuration, it suflices to show that for all x with 
|x| = 9 and y(x) / 0, e<^+i>(y(x)) = 0. However, if |y(x)| is odd, then |x| and 
T+l have different parity. Hence 7r(x, T + 1) = and (x, T + 1) = 0, and 
by construction, e^^+^^(y(x)) = 0. Moreover, for all x with |x| = 0, 

Z^^'+i^ (x, T + 1) = (x, T + 1) + e<''+i> (y (x)) = 7r(x, T + l) (mod fc) 
and for all x with |x| < 0, 

Z^^'+i^ (x, T + 1) = (x, T + 1) = ^(x, T + l) (mod fc). 

Therefore, is fine for ah (x,r + 1) with |x| < 6* + 1. 

For ah vertices x and all 9 > max(0, |x| - T), e<''>(x) = e<''+i>(x). Hence 
e(x) := lime^ooe<''>(x) and /(^+i)(x,0) := /(^'(x,0) + £(x) fc-^+i arc well- 
defined for all X € Also, /C^+i' is fine for x {0, . . . , T + 1}. 

We now note that for all vertices x £ V and for all T > |x|, we have 
that /(^)(x,0) = /(^+i)(x,0). Hence again, /(x,0) := limT_>oo /^"^Hx, 0) is 
well-defined for all x e Therefore, / is fine for F x No, which finishes the 
proof. □ 
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4 The Basic Method 



In this section, wc lay the foundations for our analysis of the maximal possible 
single- vertex discrepancy. Closely following the arguments of we will see 
that it is possible to determine the contribution of a vertex to the discrepancy 
at another one independent from all other vertices. 

For the moment, in addition to the notations given in Section [21 we also use 
the following mixed notation. By E{x, ti,t2) we denote the (possibly fractional) 
number of chips at location x after first performing ti steps with the Propp 
machine and then t2 — ti steps with the linear machine. 

We arc interested in bounding the discrepancies |/(x, — E{x,t)\ for all 
vertices x and all times t. It suffices to consider the vertex x = 0. From 

^(0,0, i) = E{0,t), 
EiO,t,t) = /(0,t), 

we obtain 

f{0,t)-E{0,t) = Y,iEiO,s + l,t}-EiO,s,t)). 

s=Q 

Let |x| denote the distance of a vertex x to 0. Now E{0, s + 1, t) — E{0, s,t) = 

ExGV EfcT"'* + next'^-i (arr(x, s))\,t-s-l)-H{\x.\,t-s)) motivates 

the definition of the influence of a Propp move (compared to a random walk 
move) from vertex x in direction A E { — 1,+!} on the discrepancy of (t time 
steps later) by 

iNF(a;, A, t) H{x + A,t - I) - H{x, t). 

In order to ultimately reduce all ARRs involved to the initial arrow settings 
ARr(-,0), wc define Si(x) := min {u > \ i < I]JLo/(^'*)} * ^ I^o- 

Hence at time Sjfx) the location x is occupied by its i-th chip (where, to be 
consistent with [3, we start counting with the 0-th chip). 

Consider the discrepancy at at time T . Then the above yields 

/(0,T)-£;(0,T) = ^ inf(|x|,nextXarr(x,0)),T-s,(x)). (3) 

xey i>o, 
s»(x)<r 

Since the inner sum of equation ([3]) will occur frequently in the remainder, 
let us define the contribution of a vertex x to be 

con(x) := ^ inf(|x|,next'(arr(x, 0)),r - .Sj(x)), 

j>0, 
s,(x)<T 

where we both suppress the initial configuration leading to the Sj(-) as well as 
the run-time T . 

We summarize the discussion so far in the following theorem. It shows that 
it suffices to examine each vertex x separately. 
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Theorem 2. The discrepancy between Propp machine and linear machine af- 
ter T time steps is the sum of the contributions COn(x) of all vertices x, i. e., 



5 Divergence of the models 

In this section, we analyze a specific initial configuration and show that the 
Propp machine may deviate from the linear machine by an arbitrarily large 
number of chips. 

We choose the initial directions of the rotors to point inwards in the direction 
of the origin. For a fixed time T at which we aim to maximize the discrepancy 
/(0,T) — E{0,T) we examine a configuration in which all vertices x with < 
|x| < T/X and A ;= are occupied by a number of chips not divisible by k 
only once. We assume that at time T — <|x| with t\^\ := [[A|x|]0an "odd chip" is 
sent in the direction of the origin. That is, the number of chips at x at time t\-^\ 
is /(x, t|x|) = 1 (modfc) while f{'x.,t) = (modfc) for all times t ^ For the 
given initial direction of the rotors such a configuration exists by Theorem [TJ 
We will prove the following theorem. 

Theorem 3. For any initial direction of the rotors and any T > 0, there is 
an even initial configuration such that the single vertex discrepancy between the 
Propp machine and linear machine after T time steps is n(\/kT). 

By Theorem[2l the discrepancy at the origin at time T of the above described 
initial configuration is 



/(o,r)-i?(o,r) = ^coN(x). 



f{0,T)-E{0,T) 



^ (if(|x|-l,t|,j 



l)-i?(|x|,t|,|)) 



x6V, 
|x|<T/A 

LT/AJ 

^ fc(fc-ir-l(i/(.T 



H{x,t,)) 




LT/AJ 



(fc-1) 



x-1 




x=l 



LT/AJ 



(fc-1) 



x~l 




(4) 



x = l 



with 



i{x, t) := k n{x — 1, t — 1) 



n(x, t) 



■"-With |yA|x|]] wc denote the smallest integer t with t > A|x| and x ^ T — t. 
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for all x,t > 1. Let us define i{x, t) = otherwise. For x,t > 2 we get 

t) = k n(x — 1, t — 1) — n(x, t) 

= kn{x - 2, t - 2) + fc(fc - l)n{x, t - 2) 

- n{x - 1, t - 1) - (A: - l)n{x + 1, t - 1) 
i{x- l,t - 1) + (fc - l)i{x + 1). 

It remains to examine the cases {x,t) £ (N x N) \ (N>2 x N>2)- Right from the 
definition, we get i(l, 1) = fc — 1. For x > 2 and t = 1 we have 

i(x,l) = kn{x — 1,0) — n{x, 1) 
= 

= i{x - 1, 0) + (fc - \)i{x + 1, 0). 

Also for a; = 1 and t > 2 we have 

i{l,t) = kn{0,t- 1) - ??.(!, t) 

= (fc - 1) (?i(0,t- 1) -?i(2,i- 1)) 
= (fc-l)(fcn(l,t-2)-n(2,t-l)) 
= i{0,t-l) + {k-l)i{2,t-l). 

Summarizing the above, we see that i(x,t) can be defined recursively as 
follows. 

i{x,0) = forallx>0, 
i{0, t) = Q for all t > 0, 

z(l,l) = fc-l, 

i{x,t) = i{x - l,t - 1) + {k - 1) i{x + l,t - 1) for {x,t) eN|i\{(l,l)}. 

This recursive view of i{x,t) reveals another interpretation of these quan- 
tities. Apart from a factor (fc — l)*^^ ^-^K, i{x,t) counts the number of lattice 
paths from the origin (0,0) to {x,t) of steps (+1,-1-1) and (—1,4-1) which do 
not cross the line a; = 0. This can be described by the well-known Ballot num- 
bers. The classical description is as follows. Suppose A and B are candidates 
for president. Let A receive a total number ofa votes and B one of b votes. Let 
a > b. Now consider the progress of counting the votes, one after the other. 
Then the probability that throughout the counting B never has more votes than 
A is {a — b + l)/(a + 1) 0. This implies that for given positive integers a,b 
with a > b, the number of lattice paths starting at the origin and consisting 
of a upsteps (4-1, +1) and b downsteps (+1,-1) such that no step ends on the 
X-axis is ("q • We are interested in the number of lattice paths starting 
at (0,a;) and consisting of {t — x)/2 upsteps (4-1, 4-1) and {t + x)/2 downsteps 
(-f 1, —1) such that no step ends below the x-axis. Therefore, 

^(x,^) = (fc"l)('-"+^)^(4) (5) 
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for x,t > 0. Note that this can be contmued to all {x,t) G 

Here is a more intuitive interpretation why the Ballot numbers come into 
play. We want the effect of having one chip at distance x and time t make 
its first move toward the origin. So we can reduce to comparing two chips, one 
at distance x — 1 and the other at distance a; + 1 at time t — 1. 

We consider the two chips on a common random data set. The data set is a 
string s of length t — 1 from {/, O} with probabilities ^, for / (in, toward 
the origin), O (out, away from the origin) respectively. The string s determines 
where the chip goes (in terms of how far away it is from the origin, which is 
what concerns us): with O the distance x goes up by one and with / it goes 
down by one unless the current distance is zero in which case it goes up by one. 

Say a string s has the property MINUS if for some initial segment the Ps 
exceed the O's by precisely x. In a "fictitious walk" , allowing distance x to go 
negative, the chip starting at x would reach —1. But consider just before this 
happens for the first time. The two chips are at 0, 2 respectively and an / is 
received. Now the two chips are both at 1. Hence whatever comes after, they 
end up at the same point. Thus the two chips have the same probability of 
MINUS and ending at zero. The times that the near chip ends at zero but the 
far chip does not are then precisely the strings s which do not have MINUS and 
for which the /'s exceed the O's by precisely x —1. And as the far chip is never 
closer than the near chip, it never happens that the far chip ends at zero but 
the near chip does not. Thus the difference in the probabilities is given by the 
Ballot problem as calculated above. 

Equations (|4]) and ^ now give 



LT/AJ 



/(o,T)-£;(o,T)= J2 (fc-i) 

(fc_l)(^) X f t 



> 



x = l 
lT/\ 

LT/AJ 

"2AP^ 



(fc-1)^- / WXxl . 

2^ o\h^^ \\M±£r ^' 



x = l 



It remains to bound the binomial coefficient. From the Stirling formula we know 
that 

for 71 > and get 

Ai\ 71"+!/^ 
3 (n - /c)"-'^+i/2 \k) 2 (71 - fc)"-'=+i/2 fc'^+i/2 ■ 

for n,k > 0. Note that equation (|5]) also holds for the generalized binomial 
coefficient which is defined for nonintegral ti, fc by using the Gamma function F. 
The following lemma shows that inequality ([6]) also holds without the ceiling 
function in the binomial coefficient. 
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f 2x \ 

Lemma 4. ) is monotonic nondecreasing in x for < y < x. 

\x + yj 

Proof. With denoting the logarithmic derivative of the gamma function, i. e. 
the Digamma function, we get 



d 



2x 

x + y 



dx 



(2*(2a;+l)-*(.T + y+l)-«'(a;-?/ + l)) ( ). 

' \x + yj 



The lemma follows by the fact that above binomial coefficient is positive for 
a:, y > and that '^{x) is monotonic nondecreasing for a; > (0, Thm. 7]). □ 

By the definitions of A and tx wc now get with Lemma |4] and equation ([8]) 



for x > 0, fc > 2. 

Using this we obtain for all fc > 2 

/(0,T,-E,0,T,> g _j(_il_.o(VEr), 

which proves Theorem [3l Using the same arguments as the following two sec- 
tions one can also prove that for all even initial configurations the single vertex 
discrepancy after T time steps is at most 0{^/kT). 



6 Convergence of the models 

The previous section showed that for very special configurations the single vertex 
discrepancy can be unbounded. In this section wc show that, on the other hand, 
many configurations have a bounded discrepancy. Wc will prove the following 
theorem. 

Theorem 5. // /(x, t) = (modfc) for all x and t such that (1 — £)A|x| < 
T — t< (l + e)A|x| with A :~ 77^27 then the discrepancy between Propp machine 
and linear machine at time T and vertex is bounded by a constant depending 
only on k and e > 0. 

Our aim is to bound equation ([3]). To that end, we further examine INF. 
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From its definition and equations ^ and (O we know for x,t > 

i{x, t) 



iNF(a;, —1, t) 
lNF(a;, +l,t) 



k n{x + 1, i — 1) — n{x, t) 
-^n{x, t) - -j^n[x + l,t-l) 



— i(a;, t) 
^ (fc- 

Tliis also shows 

lNF(a;, -1, t) + {k- 1) inf(x, 1, t) = 0. (10) 

Therefore, the absolute value of the influence |iNF(a;, A, t)| of sending one chip 
towards is (fc — 1) times larger than sending one chip in the opposite direction. 
Note that 

coN(x) = y 

i>0, 
s.(x)<T 



with 



for NEXT* (arr(x, 0)) = +1 
1 for NEXT* (arr(x, 0)) = -1. 



To bound this alternating sum, we use the following elementary fact. 

Lemma 6. Let / : X — > M be non-negative and monotone nondecreasing with 
X CR. Let A^°\ . . . , e R and to,..., t„ € X such that to < ... <tn and 
I ELa \<l for allQ<a<h<n. Then 



i=Q 



< max/(a;) 



Let X C R. We call a mapping f : X ^ M. unimodal, if there \s a. ti £ X such 
that f\x<tx as well as f\x>ti are monotone. The following lemma shows that 
lNF(a:, A, t) is unimodal. 

Lemma 7. The function i{x,t)/k^ is unimodal in t for all x and k, and it is 
maximized over all t € N at some tmax(2;) — Xx ~{- Cx with \ca:\ < 15. 

Proof Using (l+l) = (l) we get 

i{x,t + 2) i{x,t) {k-l)^xp4t) 



fc*+2 fc* {t + x + 2) {t~x + 2) A:*+2 i 



t+x 
2 
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with p^(t) := - 5fc2 + 8fc - 4)*^ - (4/fc3 - 8k^ + 8k- 4)t + (fc^ - /t2)(a,2 _ 4)^ 
Hence, the above difference is non-negative if 

t > ( V8fc3 - 4fc2 - 8fc + 4 + F (fc _ 2)2 1-2 - 2fc2 + 2fc - 2) / (fc - 2f (11) 

and non-positive otherwise. Thus we have unimodahty with i(x,t)/k* taking 
its maximum when t is the smallest integer having the same parity as x and 
satisfying equation (fTT]) . Using fc > 3, the claim follows. □ 

Armed with Lemmas |6] and [7l we can now prove Theorem [5] 

Proof of Thru. [5l By Lemmas [6] and [7] we know 

z(|x|,L(l-e)A|x|J) , z(|x|,r(l+e)A|xn) 



con(x) < 



fcL(l-e)A|x|J /,r(l+e)A|xn 

By Theorem [2] it therefore remains to show that 

/(o,r)-i?(o,T) 

^E(pfeSl^^(-' L(l-e)AxJ) 



x>0 



(fc - 1 

"fcr(i^ 



is bounded. We now show that the second summand is bounded, the first one can 
be handled analogously. For this, we choose e' such that (l+e')Aa; = [(l+£)Aa;] . 
By equation ([5]), 

|i^i£^^(-'(l+-')Ax) 

l)((i+c')A+iW2^ / (i + e')Ax 



k{l+e')Xx-l(^(^l + i^f^xx) V((l + £')A + 1) x/2 
With equation ^ we get. 



/ fc(fc-2)3 

V (1 + + 2)(2 fc - 2 fc-e') ' ' 



with 



, „ 2fc - 2 -f ke' 
p{e ) : 



(fce' -I- 2)(fc - 1) V (fe' + 2)(2fc - 2 + fc-e') 

-l)(2 + 2£') 
2fc - 2 + fce' 



ke'/2 



'(fc - l)(2 + 2£')"'+'V 
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As < p(e') < 1 for all e' > and fc > 3, we have shown that 

Sx>o fc(i+E')Ax-i (1 + s')^^) is bounded above by a constant depending on 

e' > 0. The same arguments demonstrate that X]a;>o rfe^^'^feFT *('^' L(l ~ 
e)Aa;J) can also be bounded above by a constant depending on e > 0, which 
finishes the proof. □ 



7 Number of Chips 

In the previous sections, we have shown that the discrepancy on a single vertex 
can be arbitrarily high provided we use a sufficient amount of time. However, a 
closer look at the proof also reveals that, assuming k to be constant, to obtain a 
discrepancy of D, at least exp{il{D^)) vertices have to contribute (by holding a 
number of chips that is not a multiple of at a suitable time). We now discuss 
that also exp(17(L)^)) chips are necessary to obtain a discrepancy of D. 

Theorem 8. Let k > 3 be a constant. Then all single vertex discrepancies 
arising from an even initial configuration with k chips are bounded by 0(^log k). 

Proof. Let us examine the discrepancy at the origin after T time steps, starting 
with an arbitrary even initial configuration which uses only k chips. We bound 
the contribution of each sphere Sx ■= {x G V \ |x| = x}, x G No, separately. 
Lemma [7] shows that lNF(a;, —1, •) is unimodal. By Lemma[6j this gives 

COn{Sx) := ^ COn(x) 
xes^. 

< ^ maxiNF(|x|, — 

= fc(fc -1)^-1 max %^ (12) 

By equation ([5]) and Lemma [71 which showed that i{x,t)/k^ is maximized at 
t = \x -\- Cx with \cx\ < 15, we obtain 

i(xA) i(x,Xx + Cx) 
max — — — = 

= X fc— (k 1)^^-^+! (A.. + Cx)-' ( f^xtl ) . 

Since both \cx\ and A = k/{k — 2) are bounded by absolute constants, we have 
Xx + Cx \ / f Xx 

{\+l)x+c^ / {\+l)x 



(Xx + cx)- {Xx + 1) 



[(A±1^+^^).....(A±1^ + 1)] . [(A_l^+i^).....(A_l^+l)] 



0(1) 
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Now equation © and using that k is seen as constant yields 

i{x,t) " -^-^^-^ 

max — — — — O 
t ¥ 



= O 




= 0(x-i/2(fc-l)--). (13) 

From equations and we conclude 

COn(5,) = 0(.T-i/2). (14) 

For large x, there are not enough chips such that each vertex on Sx contribute 
fully to the discrepancy. We use this to obtain a second bound for CO^^{Sx)- 
Theorem [S] implies that the sum of all contributions CON(S'a;) for all time steps t 
with T — t> ^Xx 01 T—t < I Ax is bounded by a constant Ck- Hence it suffices 
to examine times t with ^Xx < T — t < ^Xx. In every time step, each chip in 
Sx can contribute only at most 0{x^^^^k^^) by equation ([T3| . As there are k 
chips and Aa; possible time slots, the total contribution of Sx in the time interval 
^Xx <T -t < |Aa; is 

0(rai/2fc-"). (15) 
With Theorem [21 the two bounds of equations (fT4|) and (fT5|) now yield 

fiO,T) - EiO,T) <Y,CONiSx) 

x>0 

CON(S'^)+ CONiSx) 
l<i:<logfc(K) 3;>logj.(K) 

< J2 0(.T-l/2) + Cfc + 0{KX^'^k~^) 

l<x<log^(«;) 2;>Iogfc(K) 

= 0(71^) + o( VT^), 

still assuming fc to be a constant. □ 



8 Conclusion 

In this paper we showed that fc-ary trees (k > 3) do not admit a constant bound 
on the single vertex discrepancies. Nevertheless, also on these trees the Propp 
machine is a very good simulation of the random walk. For simultaneous walks 
of K chips, the discrepancies are bounded by 0(-\/log k). 

With this work showing that infinite regular trees do not have discrepancies 
bounded by a constant, but previous work of Cooper and Spencer [l| showing 
this property for higher-dimensional grids, the natural open problem arising 
from this work is to give more insight to the question of which graphs display 
the one or the other behavior, or ideally, to give a characterization of those 
graphs having constant bounds for the single vertex discrepancies. 
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